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Abstract 

Theory of the Casimir effect for a flat graphene layer interact- 
ing with a parallel flat material is presented in detail. The high- 
temperature asymptotics of a free energy in a graphene-metal system 
coincides with a Drude high-temperature asymptotics of the metal- 
metal system. High-temperature behavior in the graphene-metal sys- 
tem is expected at separations of the order of 100 nm at temperature 
T = 300K. 



1 Introduction 

The quasiparticles in graphene[l] obey a linear dispersion law u = vpk 
{vf ~ c/300 is a Fermi velocity, c is a speed of light) at energies less than 2 
eV. Graphene's 2 + 1 - dimensionality and quasi-relativistic Dirac model for 
its quasiparticles make it possible to derive the Casimir effect properties of 
graphene systems from general constraints and principles of quantum field 
theory. 

Casimir effect in graphene systems was studied in different papers [2]- 
[8]. Finite-temperature results were obtained in Refs. [1] and [5]. In the 
current paper we follow the formalism developed in Ref. [5], an alternative 
derivation of the reflection coefficients is presented, the relation to Feynman 
diagrams is discussed. First we derive the expressions for the components 
of the polarization operator of quasiparticles in a graphene layer at finite 
temperature and derive the refiection coefficients of a flat graphene layer 
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from the solutions of the boundary problems for vector potentials. The free 
energy is given then in two equivalent forms: in terms of reflection coefficients 
and closed Feynman diagrams. Finally we study the exceptional properties 
of the free energy of a flat graphene layer - parallel flat metal system at finite 
temperatures. 

We use the coordinates and z interchangeably throughout the paper. 
When needed we select the coordinate y along a direction of the wavevector 
q = (g\ q^) (longitudinal direction) and the coordinate x along a transverse 
direction. We use h = c = ks = I- 



2 Action and polarization operator 

The model is described by the following classical action (assuming graphene 
plane lying at x^ — 0) 

S^-ljd'x + j d'x^lp^ (1) 

with 

^ = (i5o - // - eAo)7o + vp^y^ {idi - eAi) + 7^ {id2 - 6^2)] - m . 

Here // is a chemical potential, m is a mass gap of quasiparticle excitations, 
1/300 is a Fermi velocity. Since there are N — A species of fermions in 

graphene, the gamma matrices are in fact 8x8, being a direct sum of four 
2x2 representations (with two copies of each of the two inequivalent ones), 
7q = —(7^'^)^ = 1. The Maxwell action is normalized in such a way that 

= 47ra = ^. (2) 

In Minkowski space the one-loop polarization operator can be expressed 
in momentum space as 

n"^"(po, P) = ^e' J tr [s{qo, ^)rS{qo - Po, q - p)7") , (3) 

where the propagator of the quasiparticles in graphene reads 

A. {qo + lx)-fo - VF<^ - m 

Siqo,ci) = jP k-o--(^^^^^.^^g^^^^2_,2q._^.- (4) 

Note that for /i = m = the pole of the propagator yields the linear 
dispersion law for quasiparticles in graphene: go = ^'fIqI- Here q = (g^g^). 
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^ = l^Qi + 7^Q'2- Vectors with tilde are rescaled by multiplying the spatial 
components with vp, i.e., = rj^p^ = {pq,vfp), rj = diag{l,VF,VF)- 

To introduce the temperature in ([3D we perform the rotation to the Mat- 
subara frequencies 

/oo 
dqo -2nT ^ , go ^ 2mT{k + 1/2), (5) 

k = — OD 

use the Feynman parametrization 

1 dx 



ah Jq {xa + (1 — x)by 

and subsequently change the variables in ([3]) in the spatial part of the loop- 
integration: q — )■ q + xp. Then we come to 



oo „i „ 

noo ^ _2e^TN J2 J "^^J 



d^c[ + {qok + /i)(gofc + /i - Po) 
(2vr)2 [{q,, + ^-xpoy-e^f 



AaTN ^ /•+°° ((q^ +pIx{1-x)-2vIp'x{1-x] 
7, — > I dx I dQQ\ — ;5 



(6) 



[{qok + iJ- - xpoY - 62]^ 

where Mq = + VpC[^ — x{l — x)vpp'^, and 

= m^ + vlq^ -x{1-x){pI-vIp^), (7) 



Qo = ^m"^ — x{l — x){pI — Vpp"^). (8' 

In equation ([HD we introduced the integration variable 0. 
In analogy we get 

- 2(gofc + /i)(gofc + -Po) 



+ = -AaTN ^ dx dQQ 

k=-oo 



iBo [{qok + /i - xpoY - 62] 

(9) 

Summation over the fermion Matsubara frequencies can be made explic- 
itly by making use of the identities 

oo 

E 



[{2mT{k +1/2) - by -Q^Y 
^ (e .ech^ (^) - 2Tta„h | ) + (0 ^ -8), (10) 
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{27riT{k + 1/2) + fx){2TriT{k + 1/2) + fi- po) 

k=—oo 



^ [{2'KiT{k + l/2)-hf -Q^f 



IQQ3T2 I ^^^h^ (©' + - l)po - epgx(l - X)) + 

+ 2rtanh(^) (e2+p^x(i-x))j +(e^-e), (ii) 

where in flTTl) we substituted 6 = po^; — /i- 

To perform the integration over G it is convenient to use the identity 
ds tanh s = sech s. Finally we arrive at the following representation for IIoo 
and Utr = n;;^: 

2NaT [\ f , ^ ^Go + fe , u ©o + &\ ^ ^\ 

ntr.oo = — j dx I ftrflo tauh In I 2 cosh 1 + (Bo ^ -©o) 1 

(12) 

where 6o = a/tti^ — x(l — x){pl — Vpp"^), b = pox — jj,, and 



/OO — , Uoj 



/tr 



-2t;2 p2x(l - x) - j9o(l - 2x)eo + 29^ 

2m'^Vp + 2x(l — x)vppl 

Po(l-2t;^)(l-2a:)-2(l-4)eo 
4T 



(14) 



Here = Pq — p^. We remind that N is the number of fermion species, 
N = 4 for graphene. Parity-odd contributions to the polarization tensor 
cancel out between different species, while the parity-even contributions add 
up. 



3 Reflection coefficients 

In this section we find reflection coefficients for transverse electric (TE) and 
transverse magnetic (TM) modes. The equations 

d^F^" + 5(;2)n^Mp = (15) 

lead to the conditions 

dzAm\z=+0 ~ dzAm\z=-0 = ^mnA'^\z=0, (16) 
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where [5] 

^ }^r^^(^ui'A{po,p)+plKB{po,p))v? (17) 

= 6^^ and 2, j = 0, 1, 2. Here A, S are scalar functions. 
Let's consider the condition 

doA^ + d^.A' + = 0. (19) 

In fact, the condition fll9l) is quite convenient for a description of transverse 
electric and transverse magnetic modes of the propagating electromagnetic 
wave. 

A nonzero Ax-, the condition dxA^ = and the conditions Ay = A^ = 
= describe the propagation of a TE electromagnetic wave (the electric 
field is parallel to the surface z = 0) since ~ A^. Here the direction 
perpendicular to the wave vector of the electromagnetic wave under consid- 
eration iS^.Py.Pz) is denoted by x. 
For the TE wave we have: 

Ax = e'Py^e'P^' + rTEe'^^^e-'P^' for ^ < (20) 

Ax = e'P^'e'PyHTE for z>0 (21) 

and 

UxnA^ = AxA{p). (22) 

Here pi = pi — Py- From the continuity of potentials at 2; = we obtain 
1 + rxE = txE- Now one substitutes (|2UI) and (PT]) into (ITB]) and uses to 
obtain: 

The conditions A^. = = 0, po^o = Py^y describe the transverse mag- 
netic (TM) wave. This choice of vector potentials describes TM wave since 
Ez ~ OzAq or Bx ~ dzAy. For Aq we have: 

Ao = e'Pyye'P^' + rAoe'P^ye-'P^' for z<0 (24) 

Ao = e^^'^^e^^'^^tAo for z > 0. (25) 

and 

Ho„A" = {A{p) - B{p)vIpI)pI/pIAo, (26) 
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where p1= Pq — v^Py. One substitutes and into (fTBl) and uses 
to obtain: 

Since i^^ ~ f^z^o the reflection coefficient for the TM mode is equal 

^ PzjAjp) - plvlBjp)) 



The reflection coefficients fl2^ and (1251) are reflection coefficients of trans- 
verse electric and transverse magnetic modes respectively. 

One can also rewrite the reflection coefficients in terms of the polarization 
tensor components 

P.noo _ pIUqo + pIu,, 

4 QED point of view 

The two conditions follow from gauge invariance: 

PoUoo{po) -pyUyo{po) = 0, 

Po^OyiPo) -Py^yyiPo) = 0, 

which yield after Wick rotation 

plUoo{ipo) = -plUyy{ipo) (30) 

and the property 

^tr{ipo) = noo(^Po) ^° - Uxxiipo)- (31) 

Py 

The reflection coefficients can be rewritten in the form: 



Po + Pl^yyiiPo) f y + Pl'^yyim 



rTM(«Po) = 1 7r2 1 ' (32) 



= ( 1 _ JW^^ (33) 



'2\pI + pI\ '2\pI + p 



6 



In the gauge Aq = the longitudinal part of the free photon propagator 
has the form 

the transverse part of the free photon propagator has the form: 



Po + Pye 

D\^Po) = ^ , (34) 



p-\^\Vpo+pi 

D^{ip,) = , (35) 



'i^Pi + Pi 

Let's choose coordinate axes in the plane of a graphene sheet so that 
Py = |p|. Lifshitz free energy [TT] has the form: 

^-^Y. |^lW(l-e--v^r<'>r?^)(l-e--V^r«ra)l, 

n=— OO 

(36) 

where Un = 27mT are Matsubara frequencies, the respective transverse mag- 
netic and transverse electric reflection coefficients from two parallel fiat sur- 
faces separated by a vacuum slit a are denoted by vtm and tte- For the 
ideal metal ttm = 1, t^te = —1- Note, however, that some exact results in 
complicated geometries [12]- [13] can be essentially different from the approx- 
imations based on Lifshitz formula for two parallel plates (see also Ref. [13] 
which considers sphere-plane system). 

From the comparison of the formula fl36l) and expressions fl32|) - fl35|) it 
follows that for graphene - ideal metal, graphene - graphene the Lifshitz 
theory takes into account the set of closed Feynman one loop diagrams re- 
sponsible for interaction between the two materials separated by a vacuum 
slit. The longitudinal and the transverse parts of the photon propagator 
enter two different sets of closed one loop diagrams for the free energy and 
multiplied in these Feynman diagrams by the longitudinal and the transverse 
components of the polarization operator Ilyy and Ila^a: respectively. Lifshitz 
type formulas result from the sum over closed Feynman one loop diagrams 
with 2 = a or 2; = 2a in photon propagators connecting the two sheets of 
graphene or a graphene sheet interacting with an ideal metal respectively. 
Inside each graphene layer the sum of RPA diagrams is taken into account 
by factors in round parentheses in f l32|) . fl33l) . The division of the free energy 
into longitudinal and transverse parts in terms of respective parts of photon 
Green's functions and polarization operator is equivalent to a division into 
TM and TE parts described by the reflection coefficients ttm and tte in the 
Lifshitz approach. 
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5 High-temperature asymptotics 



Let's assume m = 0, /i = [TO]. For |p| — one gets: 



trn ipo = - n°° zpo = = — ^ 

and reflection coefficients have the form: 

aNv1^\p\ 



rTE{ipo = 0) 
rTM{iPo = 0) 



aNvl\p\ + 12T' 

2aArT In 2 + aNp'^vj,/{2AT) 
2aNT\n2 + aNp'^v^p/{24:T) + \p\vy 



Zero Matsubara TM and TE terms yield the following high-temperature 
behavior of the free energy fl3^ in graphene - ideal metal system: 

" "16^ + 32aiV7r(ln2)a3 + ' " " ' ^^^^ 
•^°^'^-~192;^ + ---- 

Here 

- F I - ^ F I m) 

— — = Drude T-i-oo — TT-' id T-i-oo- [0^ 

loTra^ 2 
is the high-temperature asymptotics of the metal - metal system with a 
Drude model of the permittivity used [l5]-[25], which is equal to one half of 
the high-temperature asymptotics in the metal - metal system with the ideal 
boundary conditions or the plasma model of the permittivity used [2S] (see 
Ref . [27] for a review) . The zero frequency TE Matsubara term is suppressed 
by a factor aNvp and additional power of l/(Ta). 

The typical region of validity of the high-temperature asymptotics for the 
metal-metal system is AirTa ^ 1. However, due to a suppression of nonzero 
Matsubara terms by the coupling constant a the free energy in a graphene - 
metal system approaches the high-temperature asymptotics — T^(3)/(167ra^) 
much quicker and at shorter separations than in the metal - metal case (see 
next section). 

Note that in the graphene - ideal metal system the high-temperature 
asymptotics is derived from the first principles of quantum field theory. 
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6 Nonzero Matsubara terms 



It is often desirable to have an accurate analytical approximation of the exact 
result at different separations. We present such an expression for the sum of 
nonzero Matsubara terms in this section. 

To obtain an appropriate analytical expression we first note that at sep- 
arations H ^ vp one can put = in any nonzero Matsubara term. It is 
possible due to the exponential factor in the Lifshitz formula which effectively 
restrains the integration over impulse to a|p| < 1. In this case contribution 
of the type of f|,(a|p|)^ can be neglected compared to (apo)^ = i^nnaTY 
due to the smallness of the parameter vp. 

In the finite temperature sum of nonzero Matsubara terms in the Lifshitz 
formula one can use the reflection coefficients taken at zero temperature. The 
corrections due to finite temperature are suppressed for nonzero Matsubara 
terms, so we neglect them in the leading approximation. 

Under two mentioned above approximations and the condition m = fi = 
the reflection coefficients of a single graphene layer at zero temperature have 
the form: 



^T=o ^ 'KaN^ypl + p2 TcaN^/pl + 

™ TraN^pl + p2 + 8v/po + ^fP^ ~ T^aN^pl + p2 + 8|po| 



r=o _ TcaN^/pl + vj,p^ 7[aN\po\ 

'"te — ^ 



N^/pI + vIp^ + 8^ypl + p^ 7raN\po\+8,/pl + p^ 



(40) 
(41) 



Due to smallness of the reflection coefficients (both being of the order 
of a) we can take just the first term in the expansion of the logarithm in 
the Lifshitz formula. Note, however, that expansion of the TM reflection 
coefficient in a is not legitimate, as will become evident below (see ( H6|l ). 

The sum of nonzero Matsubara terms in f l36|l in this approximation in 
the TM case with r^^l^ = r^M° ( l40|) and r^^^i = 1 equals to 

A-Ftm = --5^[_ c^^i;:^^^r;V77;;^exp(-2as 
T 



n=l 

+ 00 



2vr^yHn/2 s, + 16nT/{aN) 

) 

exp{—Hn){l — gn + Hn) + {gnY exp(5fn)Ei((7?2 + Hn)^ 

n=l 

(42) 



lira? 



here Si = ^Juj'^ + |pP and g = 32Ta/{aN), Ei stands for the standard 
exponential integral function. It is convenient to reexpress the result fj42|) 
in an integral form. For this purpose one has to differentiate fj^2|) over H, 
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Figure 1: Ratio pi of the free energy to the leading high-temperature asymp- 
totics —T({3) I (IGvra^). Both graphs are evaluated for T = 300K. In graphene 
the values m = = were used. 




Figure 2: Ratio p2 of the free energy for a graphene - metal system with 
/i = m = to the ideal metal - ideal metal free energy at T = 300K. 
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assuming H as an independent parameter for the moment, calculate the sum 
over n and then integrate back over H (the integration constant is fixed as 
zero at — oo). Thus one obtains 



TaN /•+-^^ exp(t)(exp(t) + l)t^ 

AJ^TM = / dt ^ — n (43) 

Jh (exp(t) - 1) {8H + TraNt) 



The TE part of the nonzero Matsubara terms of the Lifshitz formula 
with the coefficients r^lj = ri^^^ from (14 ip and = — 1 gives the following 
contribution in the leading order in a 

T'^TiaN f+°° T^TiaN exp(-H) 

^^"1 ^-P(-2«^^) = - 16a (l-eM-HW 

(44) 

Thus, the complete result for the sum of nonzero Matsubara TM and TE 
terms in the approximation described above is given by (H3|) and (jSj): 

AT = AFtm + AJte- (45) 

Consequently, the leading vp = ^ contribution to the free energy is the sum 
— T(^(3)/(167ra^) + AJ-". It can be used for the comparison of the theory and 
experiment with 1% accuracy for all separations at T = 300K. 

From f H5]) and (jH]) one gets the energy at T = in the limit — )■ 0: 

where the non analyticity in a comes from the TM mode. 

FigH] clearly demonstrates that the free energy of a graphene - metal 
system approaches the high-temperature asymptotics — T(^(3)/(167ra^) much 
quicker and at shorter separations than the two metals' system. Such an ap- 
proach to the high-temperature behavior is related to the fact that nonzero 
Matsubara terms are of the order of the coupling constant a and thus very 
small in comparison with respective terms for metals. The zero frequency TM 
Matsubara term acquires the value of the sum of nonzero Matsubara terms in 
the free energy at separations a ~ lOOnm at T = 300K (see Figj2]) and dom- 
inates in the free energy at larger separations. Thus the high-temperature 
behavior in graphene - metal systems should be observed at separations of 
the order of 100 nm at T = 300K (the same effect in metal - metal systems 
takes place at separations of the order of several micrometers at T = 300K). 
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7 Conclusions 



The behavior of the free energy of the graphene - metal system is studied on 
the basis of the field theoretic model. The components of the polarization 
operator of 2 + 1 quasiparticles in a graphene layer are evaluated at finite 
temperature. The TM and TE reflection coefficients are derived from the 
solutions of the boundary problems for vector potentials. 

In the high-temperature limit the asymptotics of the free energy coincides 
with the Drude model asymptotics for two metals' system. 

The crossover to the high-temperature behavior in a graphene-metal sys- 
tem takes place at separations a of the order of 100 nm at T = 300K. This 
is the reason why the systems with graphene are very promising for the ex- 
perimental studies of the finite-temperature Casimir effect. 
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